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In this paper, we investigate a high-dimensional quantum state transfer protocol. An arbitrary 
unknown high-dimensional state can be transferred with high fidelity between two remote registers 
through a XX coupling spin chain of arbitrary length. The evolution of the state transfer is de- 
termined by the natural dynamics of the chain without external modulation and coupling strength 
engineering. As a consequence, entanglement distribution with high efficiency can be achieved. Also 
the strong field and high spin quantum number can counteract partly the effect of finite temperature 
to ensure high fidelity of the protocol when the quantum data bus is in the thermal equilibrium 
state under an external magnetic field. 
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I. INTRODUCTION 

Quantum state transfer between two distant parties is 
an important task for quantum information processing 
(QIP). The high fidelity state transfer relies on physi- 
cal systems which can serve as quantum data buses to 
connect remote parties. There are many experimental 
realization of the data buses, such as phonons in ion 
traps 0, 3 , electrons in semiconductors Q , flux qubits 
in superconductors d-jg] and photons in optics |7H10|. 
Recently, the solid-state spin chain system is becoming 
one of the most promising candidates for QIP due to its 
lon g de coherence time and ability to manipulate, transfer 
[Tl| - [l7t . The dynamics of such spin chain is determined 
by the evolution under a suitable Hamiltonian, e.g., the 
Heisenberg or XY Hamiltonian. 

The first quantum state transfer (QST) protocol was 
proposed by Bose in which the spin state can be efficiently 
transferred through a spin chain via natural evolution 
fl8j |. In the proposed protocol, two qubits are located 
at the two ends of the spin chain and the state of the 
encoded qubit at one end will be transferred to the other 
end after a specific mount of time without any opera- 
tion on the chain. In the past decades, there are many 
QIP protocols based on the spin chain system s, espe - 
cially the perfect quantum state transfer (PQST) [19l - |28j . 
Later, Christandl et al. [29, 30] generalized the PQST of 
the spin chain to the spin network with arbitrary length 
based on the Cartesian product method of graph theory. 
Yao et al. [3l|, (32[ proposed a high fidelity QST through 
an infinite temperature (unpolarized) quantum data bus. 

Entanglement and quantum parallelism in quantum 
computation provide us the potential of merits superior 
to our conventional classical methods [33]. Moreover, it is 
difficult to raise the number of qubits coupled experimen- 
tally [34| and the high-dimensional systems can be cou- 
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pled to a given dimensionality of the Hilbert space more 
efficiently using fewer systems than the two-dimensional 
systems. Hence, the high-dimensional systems can be 
employed as qudits to encode quantum information in- 
stead of qubits. The extensions of various protocols 
of quantum computation and communication from two- 
dimensional systems to high-dimensional systems have 
been proposed, such as quantum cloning [35|, |36|. quan- 
tum cryptography [37| , quantum teleportation [38|, [39j , 
quantum key distribution [40j and even implemented in 
experiment 41]. Additionally, the application of high- 
dimensional systems will enhance and deepen our un- 
derstanding of quantum computation and communica- 
tion. The high-dimensional PQST through the spin chain 
has been shown to be possible when the distances be- 
tween the two end qudits are 1,2,4 [42J and the three- 
dimensional case was realized as the distance is 3 under 
the bilinear-biquadratic Hamiltonian [43]. 

In this paper, we propose an efficient protocol to 
achieve an arbitrary unknown high-dimensional QST 
through a XX coupling spin chain. The high S spins 
are used to act as a quantum data bus and quantum reg- 
isters. The low-lying level states of registers are used to 
encode quantum information as qudits. The data bus is 
initialized to a ferromagnetic order with the spins align- 
ing in a parallel way. The states of the spin chain are 
mapped onto a set of bosons after Holstein-PrimakofF 
transformation. Under the limit that the dimension of 
the sent state is much smaller than spin quantum num- 
ber, the spin- wave interaction can be neglected to express 
the Hamiltonian of the chain in terms of free bosons and 
diagonalize the Hamiltonian of the data bus after an or- 
thogonal transformation. Ensure that the register-bus 
coupling strength is much weaker than that between the 
data bus spins, the two registers are resonantly coupled 
to one of the collective eigenmodes of the data bus and 
other off-resonant coupling can be ignored. Then a swap 
gate between the two registers can be achieved at the op- 
timal time. The proposed protocol requires neither ex- 
ternal modulation of the Hamiltonian evolution nor spin 



chain coupling engineering. Consequently, entanglement 
distribution with high efficiency can be achieved by using 
this scheme. We numerically simulate the average fidelity 
of the data bus in the thermal equilibrium state under an 
external magnetic field and find that the average fidelity 
decreases with temperature and increases with the held 
and spin quantum number. 

The paper is organized as follows. In Section II, we 
give the XX model Hamiltonian and treat it. In Section 
III, we show the high-dimensional QST with high fidelity 
and entanglement distribution with high efficiency. In 
Section IV, the case where the quantum data bus is in 
the thermal equilibrium state is investigated. And the 
last section is our summary. 



II. MODEL AND CALCULATION 
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FIG. 1: (Color online) (a) Two distant quantum registers 
are intermediated by a quantum data bus with the coupling 
strength g, which is much weaker than that between the data 
bus spins J. By using the low-lying level states of the sender, 
ranging from the |0) state to the \d — 1) state, to encode as a 
qudit, we achieve a high-dimensional QST with high fidelity. 
(b) Entanglement distribution through a spin chain. The sites 
b and c are coupled to the quantum data bus, while the site a 
not. a and b are prepared in a maximally entangled state. Un- 
der the natural evolution a and c gain entanglement from the 
maximally entangled state with high efficiency at the optimal 
time. 



A spin chain consists of N sites, where each site com- 
prises a single 5 spin in the state \m) (m — S,S — 
1,...,— 5 + 1,-5). The spin chain is used to act as a 
quantum data bus. The data bus is placed between two 
additional spins, which are considered as two quantum 
registers, denoted as the sender s and the receiver r. The 
Hamiltonian of the system consisting of spins coupled to 



their nearest neighbors on a finite lattice of site N + 2 
as shown in figure HJa) includes two terms (We have set 
H = 1 in the present paper) : 



H = H 



xx 



H 



i\/. 



(1) 



where Hxx = Hb + Hj. Hb is the XX coupling of the 
data bus as 



H B - -J} XSfSj + s i £/)> 

(id) 



(2) 



where J > is the coupling strength between the data 
bus spins, (i,j) means that only the nearest neighbor 
(NN) coupling is considered. Sf = Sf ± iSf and S*f 
[y = x,y,z) is the v component of the spin operator Si 
at the site i. The interaction Hamiltonian between the 
two registers and data bus has the form 



Hj = -g(S+Si + S+S N + H.C.) 



(3) 



where g > is the register-bus coupling strength and 
H.C. denotes the complex conjugate. Hm is the Zecman 
term under an external magnetic field h given by 



/V 



Hm — —h(S* 



£s?). 



(4) 



After the Holstein-Primakoff (HP) transformation [44 [ 






25 - alaiQi, 



a\ai, 



(5) 



with 



,ti 



Sij and [c4,al] = [o,-, Oj] = 0, the two 



Hamiltonians Hxx and Hm are expressed in terms of 
the bosonic operators. 

Initially, the data bus and receiver has a simple ferro- 
magnetic order with the spins aligning in a parallel way 
denoted by |0)fj» r with |0)fj = |0)i...|0)^. The low- 
lying level states of the sender, ranging from the |0) state 
to the \d — 1) state, are employed to encode information 
as a qudit. The total boson number, 



N,, 



N s +N r + 



N 



(6) 



is conserved. The dimension of the Hilbert space H 
associated with the spin-S 1 chain of length N + 2 is 
(25 + l) N+2 . The dimension of the state that should 
be sent is d and the conservation law ensures that the d- 
dimensional state transfer dynamics is determined by the 
evolution in the d Ar+2 -dimcnsional subspace S spanned 
by the basis vectors \n s ni...n^in r ) (n, = 0, ...,d— 1). We 
assume that the dimension d of the sent state is much 
smaller than 25, as d -C 25, which gives {a\ai) -C 25. 
The HP transformation is simplified to I45| 



S+ = ai V2S,S- 



25. 



(7) 



Substituting Eq. into Eqs. ©, © and gj), we 
have 

H B = -2SJJ2(aUj + Oi<»5)> ( 8 ) 

Fr = -2%(a|ai + aja^ + #.C.) (9) 



and 



N 



H M = -{N + 2)hS + h(a\a s + a\a r + ^a\ai). (10) 

8=1 

Here the spin- wave interaction of the system is neglected 
and the three Hamiltonians are reexpressed in terms of 
free bosons. 

As gj J <C 1, Hb works as a collective Hamiltonian and 
Hi works as a perturbation one, the diagonalization of 
the Hamiltonian Hb of Eq. J5J occurs after the following 
orthogonal transformation[3l|, |32|, |46[ 
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A ^ ~~ N + 1 

fc=l 



'k- 



(11) 



with k = 1,...,N and A = d j" 1 ■> * ne Hamiltonian Hb 
becomes 



A' 



Hb = ^ ekbjbk, 



(12) 



fc=i 



where Sk = — 4S'Jcos(-^ L j-). Combining Eqs. ^ and 
(fTT|) . it gives that 



JY 



Hi=j2 tk( - asb i + (-i) fe_1 «^i + h - c -\ 



(13) 



fe=i 



with t fc = -^ sin J$i 



III. QUANTUM STATE TRANSFER AND 

ENTANGLEMENT DISTRIBUTION 



In Heisenberg picture, the evolution of the opera- 
tor is determined by 0(t) = U'(t)O(0)U(t), where 
U(t) — e" lHT is the evolution operator. The operators 
A and B obey the rule that e - aA Be aA = B- a[A, B] + 
%r[A, \A,B\\ + .... By applying the relations we get 



u! ff aiU eff = ai + -(ai + (-ir- 1 ai) 



[- 1 + cos(v / 2i k t)] + ibl 



sin 



{V2t K 



V2 



(15) 



At the optimal time r = tq = it/y2t K: we have 

Ul ff a\U eff \ T0 = (-1T4. (16) 



Similarity, 



Ul ff atU eff \ T0 = (~lTal 



(17) 



Eqs. (fT5|) and (jTTJ) show that the creation operator re- 
ferring to the sender (receiver) at r = becomes that 
referring to the receiver (sender) at the optimal time tq. 
Actually, it is a swap gate between the two registers with 
an additional phase (— 1) K . The additional phase is in- 
dependent of the sent state and only determined by the 
length of the data bus. Thus the receiver will obtain an 
arbitrary unknown state that is sent at initial time after 
a phase gate operation given by 



pd 



/ 1 

(-1) K 
(-1) 2k 



\ 



\0 







(_1)(<*-1)k i 



Af+l 



(18) 

For simplicity, we take the 4-dimensional state trans- 
fer as an example. The sent initial state is \<p) s — 

£j=o«>> = El=o^(4r\0\ ^d El=oM = 1- 
The initial state of the whole spin chain, including the 
sender, receiver and data bus, shall be 



A. Quantum state transfer 

In this section, the high-dimensional QST with high 
fidelity through a spin chain in the absence of an external 
magnetic field is proposed. If the spin number in the 
data bus is odd, there is a zero energy bosonic eigenmode 
corresponding to k = k = (N + l)/2. By maintaining 
g/J <C 1 to ensure t K <C \s K — e K ±i|, the off-resonant 
coupling to other bosonic eigenmodes can be neglected 
and the two end re gist ers are resonantly coupled to the 
nth eigenmode [3l|, [32| with t K = —2Sg/A. The state 
transfer dynamics is drived by the evolution under the 
effective Hamiltonian 

H eff = t K {a% + (-lr^aX + H.C.). (14) 



|^(0)) 



& l4) 



10)] S |0)^|0) 



(19) 



At the optimal time To, the creation operator referring 
to the sender becomes that referring to the receiver, we 
obtain the final state 

I^(to)) 



io).io>^E-^f(-ir(4)"i 

|0).|0)g^E(-l)' Ms a M |/i>]r. 

ju=0 



(20) 



After a phase gate operation 



P, 



10 

(-1) K 

v - I (-1) 2K 

(-1) 3k 



(21) 



JV+1 



we have P^ +1 |^(r)) = |0) s |0)i ms |<^) r . 

In two-dimensional Hilbert space, the set of the pure 
states forms a complex projective space CP 1 , one can use 
the two-dimensional Bloch sphere to measure these pure 
states, since in this case U(2)/U(l) ~ 50(3)/5O(2) ~ 
52- The set of the pure states of a d-dimensional Hilbert 
space constructs a complex projective space <CP d ~ l , 
where the natural uniform measure exists. To measure 
a random pure state on this 2(d — l)-dimensional man- 
ifold, a vector of a random unitary matrix distributed 
over the invariant (Harr) measure on U(d) is taken. A 
rf-dimensional normalized pure state \<f>) can be measured 
by the Hurwitz parametrization with d — \ polar angles 
X P and d — 1 azimuthal angles 9 p as [47| 



(cos 6 d -i, sin d _i cos Od^e 1 ™- 1 , sin d -i sin 9 d - 2 



d-l 



cos0 d _ 3 e lXd - 2 ,...,n sm ^ elXl )' 



(22) 



i=l 



where < P < ir/2 and < \p < 27r with p — 
1,2, ...,d — 1. In analogy to the volume element on the 
two-dimensional Bloch sphere, that on the generalized 
Bloch sphere of Eq. (f2"2l in the complex space CP d_1 
can be written as 



d-l 



dV =Y[ cos 9 p {sin Op) 2 ?' 1 d6 p dxp 
P =i 



(23) 



the total volume of the 2(d — l)-dimensional manifold of 
the pure states is V d = TT d ^ 1 /(d — 1). The average of an 
observable O over the whole manifold of the pure states 
is (O) = fcp^ OdV/V d . 

The sent state is \(p) s — X)u=o a ^l/ i ) ano - the hiitial 
state of the entire spin chain is 



d-l 



hHo)) = 5>» a |o>®» f 



(24) 



The state at time r shall be 



IVK 7 ")) = ^2 c (n B n 1 ...n N n r ,T)\n s n 1 ...n N n r ) (25) 

(n s ni ...riNKr) 



and 



d-l 
(n s n 1 ...n N n r ,r) ~ / y a p.J(n,m ...n wn,V 



(26) 



wh ere f^ ni „ MNnr) = (n a » 1 ...n^n r |e- <ff, -|M>.|0)^|0>r. 

The state of the receiver at r is denoted by p r (r) since it 



is generally a mixed state. It can be obtained by tracing 
out the other sites 

Pt{t) 

= tr f (|V(r)>^(r)|) 

d-l 

= Y. P(nr<,T)\n r )(n' r \ (27) 

n r ,n' r — 

with 

P{n r n' r ,T) = 2.^1 C (n'^n'{...n'^n r ,T) c ( n ' : ln'{,...,n'l i n' T .,T)- 

(28) 
The fidelity between the sent initial state l^)^ and 
the received final state Pr(r) is defined by F(t) = 
s (ip\p r (T)\ip) s , which turns out to be 



F(r) 



d-l 



'A, 



/ 4 a ^ !0L (x^^r 



(29) 



The average of the fidelity over the complex projective 
space CP d_1 is 



(F(t)) 



1 



F(r)dV. 



(30) 



CP d 



Figure [5] shows the numerical results of the average fi- 
delity at the optimal time To. The average fidelity as a 
function of g/J is depicted in figure EJa). (F(t )) de- 
creases with either g/J increasing or spin quantum num- 
ber S decreasing when g/J > 1. As g/J — > 0, it reaches 
its maximum value, which increases with spin quantum 
number. Noted that (F(to)) could reach high values at 
some g/J values which are not much smaller than unity 
for high spin quantum number. The two cases where 
N = 3, d = 3 and N = 3, d = 4 for three g/J values 
are depicted in figures Efb) and [2jc) , respectively. It is 
seen that (F(tq)) increases with spin quantum number 
S. As g/J « 1 and rf « 25, (F(t )) shall tend to 1, 
for example, (F(tq)) = 0.9987 in the case N = 3, d = 3, 
g/J = 0.1, 5 = 10 and (F(r )) = 0.9946 in the case 
N = 3, d = 4, g/J = 0.1, 5 = 10. 



B. Entanglement distribution 

Entanglement has been considered as an important 
resource in QIP [33. And entanglement distribution 
among involved parties is also required in quantum tele- 
portation [48| . quantum dense coding [49| and quantum 
key distribution [50]. The proposed model can also be 
generalized to complete entanglement distribution be- 
tween two remote parties. 

We assume that the sites b and c are coupled to the 
quantum data bus but site a not as seen in figure 1(b). 
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FIG. 2: (Color online)The average fidelity at the optimal time 
ro. (a) The average fidelity as a function of g/J at N — 3 
and d = 3. The average fidelity as a function of spin quantum 
number at (b) AT = 3, d = 3 and (c) N = 3, d = 4. 



At initial time, a and & are prepared in a maximally en- 
tangled state 



_. d-l 



Vd 



b- 



(31) 



n=o 



The purpose is to transfer this entanglement by the nat- 
ural dynamics of the data bus to achieve a maximally 
entangled state between a and c. 

The state of the composite system can be described as 



i^ )) = ^Ei^-[%Ti ) b i°)^i ^ 



\fd r , 



(32) 



After evolution, the finial state at the time tq is 



iv>(to)> = io) b |o)^[-±=5>in M ) a | M ) c ] 



(33) 



^=o 



The maximally entangled state between a and c has been 
achieved. The logarithmic negativity [5l], [52J is used to 
measure the entanglement of the state p\2, which is given 
by 



LE l2 = log 2 || p x \ 



■i\ 



(34) 



Here T\ denotes the partial transpose of the density ma- 
trix p\i with respective to the subsystem 1 as 



(«i, J2|p^2 l fc i: ? 2> = (k\,j<z\pvz\hM) 



(35) 



The trace norm of operator O is || O ||= tryOUJ = 
1 + 2 1 J^i A«|, Ai is the negativity eigenvalue of p\\. For 
a pure maximally entangled state of Eq. (|3"Tj) , the loga- 
rithmic negativity is yielded that LE a b = log 2 d from Eq. 
(f3~4"| . The efficiency of entanglement distribution [42| is 
defined by 



E = 



LE ac 
LE ab ' 



(36) 



which is used to measure the entanglement between a and 
c gained from the maximally entangled state between a 
and b under the spin chain evolution. 

Figure 3 demonstrates the efficiency of entanglement 
transition through a data bus of length N = 3 as a func- 
tion of spin quantum number at the optimal time To. 
It is obvious that the efficiency increases with spin quan- 
tum number and approaches to unity when spin quantum 
number is high enough. For example, E(tq) = 0.9984 at 
d = 3, S = 10 and E(t ) = 0.9948 at d = 4, S = 10. 



1.00 




FIG. 3: (Color online)The efficiency of entanglement transi- 
tion through a data bus of length N = 3 when g/J = 0.1 at 
the optimal time tq. 



IV. QUANTUM DATA BUS IN THE THERMAL 
EQUILIBRIUM STATE 

In the following, we will investigate the effects of tem- 
perature on the quantum data bus which is in the thermal 
equilibrium state in the presence of an external magnetic 
field ft. 

In the case where ft ^ 0, the evolution operator can be 
described as 



e -iHr _ -iHxxr -iH M r 



(37) 



since [Hxx,Hm] = 0. The finial state at the optimal 
time To under the evolution operator of Eq. (|37]l starting 
from the initial state of Eq. (|24[) is 



d-X 



Mr)> = |0).|0>g£5> 



■1)"» 



-ihp,T 



• IaO, (38) 




0.70 
0.0 



FIG. 4: (Color online)The average fidelity as a function of 
temperature when N = 1, d — 3 in the cases S = 3 and 5 at 
the optimal time tq. 



m=o 



where the overall phase induced by — (N + 2)hS has been 
neglected. An additional phase e~ lhfJ,T ° is introduced by 
the field. When the quantum data bus is in the thermal 
equilibrium state, the density matrix satisfies the Boltz- 
mann distribution as 



PB 



1 

— < 
Z 



iE- 



*!&><& 



(39) 



where \<j>i) and Ei are the eigenvectors and eigenvalues of 
the data bus, respectively. Z = tr(e~ H ^ T ) is the partial 
function and T represents the temperature. The Boltz- 
mann constant ks has been set to be unity. As ft- = 0, 
the state \0)f us with zero energy is not the ground state 
of the data bus. As ft ^ 0, a magnetic field is added into 
Eq. P^|) and it becomes 



Figure 0] shows that the average fidelity decreases with 
temperature and increases with the field. With the tem- 
perature increasing from zero, the number of the excited 
bosons increases in the data bus, which could not ensure 
the condition (ajaj) <C 2S, to lower the average fidelity. 
On the other hand, the field depresses the excitation of 
the bosons to reduce the effect of temperature and raise 
the average fidelity. A feature of figure [4] is that all red 
curves with S = 5 vary not so steeply as the correspond- 
ing black ones with S = 3, which means that the high 
spin quantum number can counteract partly the effect of 
temperature. Therefore, at finite temperature the strong 
field and high spin quantum number can ensure high fi- 
delity of our protocol. 



A 



H B = ^2{s k 



h)b\b k 



hNS. 



(40) 



fc=i 



JV+l 



hNS. Hence, |0>®£ 
At zero tempera- 



The eigenvalue of the state |0)f^ s is 
is the ground state as ft > 45 J cos 

ture the data bus is in |0)^ with no excited boson, and 
the data bus spins are in a ferromagnetic order. With 
the increase of temperature, the number of the excited 
bosons becomes larger. 

The initial state of the whole spin chain is 



d-x 



P (0)= ]T apa* tl ,\ l i) B \0) r p B M 



(41) 



when the data bus is in the thermal equilibrium state. 
The state of the reduced density matrix of the receiver 
at time r is p r (r) = tr? [e~ iHr p(Q)e iHT }. Figure 0] de- 
picts the average fidelity as a function of temperature at 
the optimal time To when the data bus is in the thermal 
equilibrium state for several fields. 



V. SUMMARY 

In this paper, we have demonstrated that an arbitrary 
unknown high-dimensional state can be transferred be- 
tween two parties at the optimal time with high fidelity 
by arbitrary length. Consequently, the entanglement dis- 
tribution between two remote parties can be realized with 
high efficiency. This protocol does not require external 
modulation and coupling engineering. Its a direct ap- 
plication is to communicate between two remote regis- 
ters in a quantum computer using high-dimensional sys- 
tems, which can construct a bigger Hilbert space to pro- 
cess much more information than two-dimensional ones. 
When the quantum data bus is in the thermal equilib- 
rium state under an external magnetic field, the average 
fidelity decreases with temperature and increases with ei- 
ther the field or spin quantum number. Thus both the 
strong field and high spin quantum number can ensure 
high fidelity of our protocol. This helps the protocol to 
be able to work in the thermal environment. 
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